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(57) Abstract 

The invention relates to a magnetic resonance (MR) imaging 
method for fonning images of a part of a human or animal body 
which is arFBiiged in a steady magnetic field, an KIR image being 
made of the part which is directed transversely of a major axis of 
an ellipsoid which approximates the part- The method is used, for 
example for functional examination of a heart of a human body 
during which MR images are made of the left ventricle of the 
heart. Because the heart in every body may be oriented differently 
relative to a feet-to-head axis and relative to an axis transversely 
of the fect-to-head axis, the major axis of the ellipsoid should 
be determined so as to enable an MR image to be formed of an 
imaging plane which extends transversely of the major axis and 
contains a minor axis of the ellipsoid. To this end, two MR images 
are made of a first imaging plane and a second imaging plane 
of the heart, respectively, the first imaging plane and the second 
imaging plane being chosen in such a manner that they extend 
transversely of the feet-to-head aitis of the body. Subsequently, 
a first reference and a second reference are determined in the 
first and the second MR image, respectively, said references being 
characterized by one and the same property in the MR images. 
The estimate of the major axis of the ellipsoid is then determined 
by a connecting line through two points which correspond to the 
first reference and the second reference, respectively. 
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Method of and device for imaging an object by means of magnetic resonance. 



The invention relates to a magnetic resonance (MR) imaging method for 
imaging a part of a human or animal body arranged in a steady magnetic field, which method 
involves the use of imaging pulse sequences for the measurement of the MR signals for 
reconstructing the MR images, approximation of the part of the body by an ellipsoid. 
5 imaging of an imaging plane of the part which contains a minor axis of the ellipsoid, which 
minor axis is derived from a major axis of the ellipsoid, and which method includes the 
following steps for estimating the major axis: 

a) formation of a first MR image of a first imaging plane and a second MR image 
of a second imaging plane of the part of the body, the first and the second imaging plane 

10 being chosen so as to extend transversely ofa common axis of the body, 

b) determination of a first reference from the first MR image and a second 
reference from the second MR image, and 

c) determination of a first estimate of the major axis of the ellipsoid from a 
connecting line through two points which correspond to the ffa-st reference and the second 

15 reference, respectively. 

The invention also relates to an MR device for carrying out such a method. An imaging 
plane is to be understood to mean herein a slice in the body of which an MR image is made. 

A method of this kind is known from United States patent US 5,107.838. 
The known method is used to form MR images for functional analysis of a part of the body. 
20 for example the heart of a human or animal body. It is thus used, for example to form 
perfusion MR images of the heart or MR images of the coronary artery system around the 
heart. According to the known method, the ellipsoid is used as a model description of. for 
example the left ventricle of the hean. Furthermore, for the formation of perfusion MR 
images or MR images of a coronary artery of the heart it is important that an imaging plane 
25 of the heart is reproduced which contains the minor axis of the ellipsoid. Because in every 
body the heart may be oriented differently with respect to a common axis of the body, for 
example a head to feet axis, the major axis and the minor axis of the elUpsoid must be 
determined relative to an axis transversely of the head-feet axis so as to form the MR image 
of the imagir^ plane containing the minor axis of the ellipsoid. The known method utilizes 
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the two MR images so as to determine an orientation of the major axis, for which an 
operator determines, by way of a visual estimate, the first reference from the first MR image 
and the second reference from the second MR image. For the first reference, for example, a 
point in the first MR image is selected which corresponds to a point at an edge of a valve of 
5 the heart. For the second reference, for example, a point in the second MR image is chosen 
which corresponds to a point near a tip of the left ventricle. The major axis is subsequently 
determined by the connecting line between said points. The minor axis is subsequently 
determined from the major axis by way of an empirically determined ratio of the major axis 
to the minor axis of a human heart. It is a drawback of the known method that the estimate 
10 of the major axis is inaccurate and yields an MR image of an imaging plane of the body 
which deviates from an imaging plane of the body which contains the actual minor axis. 

It is an object of the method according to the invention to reduce the 
deviation between the imaging plane reproduced by the MR image and an imaging plane 
containing the acmal position of the minor axis. To achieve this, the method according to the 
15 invention is characterized in that the first reference and the second reference are determined 
by one and the same property. The invention is based on the recognition of the fact that the 
use of the same property for determining the first and the second reference enables the 
operator to make a better visual estimate; the use of geometric properties of the ellipsoid and 
the first and second reference offers a more accurate estimate of the major axis of the 
20 ellipsoid, yielding in a more accurate estimate of the minor axis so that the deviation between 
the imaging plane reproduced by the MR image and the imaging plane containing the actual 
position of the minor axis is reduced. Such a property for determining the first and second 
references is. for example, a center of a circular region in the MR image. Such a circular 
region may represent, for example, a cross-section of the left ventricle of the heart. 
25 A further advantage of the use of the same property in the first and the 

second MR image consists in that the method can be simply implemented in an algorithm for 
determining said property, thus enabling automatic recognition of the property in the MR 
images so that a visual estimate by the operator might be dispensed with and a reproducible 
estimate of the major axis and the minor axis is obtained. 
30 A special version of the method according to the invention is 

characterized in that it also includes steps for 

formation of a thud MR image of a third imaging plane and a fourth MR image 
of a fourth imaging plane of the part of the body, the third and the fourth imaging plane 
being chosen so as to extend transversely of the first estimate made for the major axis, and a 
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Step in which a second estinaate is made for the major axis by repeating the steps b) and c) 
while replacing the first MR image by the third MR image and the second MR image by the 
fourth MR image.. This step yields a second estimate of the major axis of the ellipsoid which 
is. a more accurate approximation than the first estimate of the major axis, thus reducing the 
5 deviation between the imaging plane reproduced by the MR image and the imaging plane 
containing the acmal position of the minor axis. 

A funher version of the method according to the invention is 
characterized in that said property is determined by a center in a circular region of said MR 
images of the imaging planes of the part of the body. It has been found that circular regions 
10 in the MR images represent cross-sections of, for example the left ventricle of the heart and 
that the estimate of the major axis containing the centers of the circular region of a first and 
a second imaging plane of the heart constitutes a suitable approximation of a major axis of 
the ellipsoid. It is a further advantage that algorithms for determining centers of circles are 
reasonably firm, which means that said algorithms yield substantially the same results. 
15 A farther version of the method according to the invention is 

characterized in that in order to make a next estimate of the major axis which replaces a 
previous estimate of the major axis, the method also includes the following steps: 

determination of an orientation line through a first point which corresponds to 
the first reference, substantially perpendicularly to the imaging planes, and 
20 - determination of the next estimate of the major axis from a combination of a 

predetermined ratio of the minor axis to the major axis of the ellipsoid and an angle between 
the previous estimate of the major axis and the orientation line. The use of anatomic 
propenies of the heart, for example the approximation of the left ventricle of the heart by 
means of an ellipsoid of revolution having the predetermined geometric ratio of the minor 
25 axis and the major axis, offers a more accurate estimate of the major axis. 

A further version of the method according to the invention is 
characterized in that the imaging pulse sequences include an Echo Planar Imaging pulse 
sequence. The Echo Planar Imaging pulse sequence is known from EP-A-604441 and 
includes a fast MR imaging pulse sequence for reconstructing an MR image of an imaging 
30 plane of the pan of the body. 

A further version of the method according to the invention is 
characterized in that in order to enhance the contrast between the parts representing blood 
and the parts representing tissue in said MR images, the method includes a step in which a 
pre-pulse sequence is generated prior to the generating of the MR imaging pulse sequence. 
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Such a pre-pulse sequence enhances the contrast between, for example parts in the MR image 
which represent blood present in myocardium and parts in the MR image which contain the 
myocardium. A T^Y pre-pulse sequence is an example of such a pre-pulse sequence. The 
TjY pre-pulse sequence is known from the article "Myocardial Suppression In Vivo by Spin 
5 Locking with Composite Pulses", published by W.T. Dixon et aL in Magnetic Resonance in 
Medicine, No. 36. pp. 90-94, 1996. An advantage of the enhanced contrast in the MR image 
resides in the fact that it enables more accurate estimation of the references from the MR 
images. 



10 characterized in that it includes a step in which a contrast medium is introduced into a blood 
stream through the part to be imaged. An example of such a contrast medium contains 
gadolinium-diethylenetetramine pentaaceticacetate (Gd-DPTA). This known paramagnetic 
contrast medium reduces the longitudinal magnetization relaxation time Tj as well as the 
transverse magnetization relaxation time To, depending on the concentration used. 
15 The invention also relates to a device for performing such a method. 

A device according to the invention for magnetic resonance imaging of a part of a body of a 
human or animal which is to be arranged in the device, which device also includes means for 
sustaining the static magnetic field, means for generating gradients, means for radiating RF 
pulses towards a body placed in the static magnetic field, control means for steering the 
20 generation of the gradients and the RF pulses, means for receiving and sampling MR signals 
generated by sequences of RF pulses and switched magnetic gradient pulses, means for 
reconstmcting a first MR image and second MR image plane, respectively of a part of the 
body to be arranged in the device, which image planes are transversely of a conmion axis of 
the body arranged in the device, means for determining a first reference from the first MR 
25 image of the first image plane, means for determining a second reference from the second 
MR image of the second image plane, means for determining a first estimate of a major axis 
of an ellipsoid which approximates the part of the body to be arranged in the device, said 
first estimate being determined by a connecting line between two points, said two poings 
corresponding to the first and second reference, respectively, and means for forming an MR 
30 image of an imaging plane of the part which contains a minor axis of the ellipsoid. 
The device is characterized in, in accordance with the invention, in that the means for 
determining the first reference and the means for determining the second reference are 
arranged to determine one and the same property from the ilrsi and the second MR image, 
respectively. 



A further version of the method according to the invention is 
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The above and other, more detailed aspects of the invention will be described in detail 

hereinafter, by way of example, with reference to the drawing. 

Therein: 

Fig. 1 shows a magnetic resonance imaging device, 
5 Fig. 2 shows an EPI pulse sequence. 

Fig. 3 shows a first cross-section of the imaging planes and the ellipsoid. 
Fig. 4 shows a first contour diagram of the second class inaccuracy as a 
function of the ellipsoid ratio and the angulation, 

Fig. 5 shows a second cross-section of the imaging planes and the 

10 ellipsoid. 

Fig, 6 shows an error diagram of a first-class inaccuracy, and 
Fig. 7 shows a second contour diagram of the second-class inaccuracy as 
a function of the ellipsoid ratio and the angulation. 

Fig. 1 shows a magnetic resonance imaging device which includes a first 
15 magnet system 2 for generating a steady magnetic field, and also various gradient coils 3 for 
generating additional magnetic fields which are superposed on the steady magnetic field and 
cause a gradient in the steady magnetic field in three respective orthogonal directions of a co- 
ordinate system X, Y. Z. Generally speaking, a gradient in the first direction is referred to 
as a read-out gradient, a gradient in the second direction as a phase encoding gradient, and a 
20 gradient in the third direaion as a selection gradient. The Z direction of the co-ordinate 

system shown corresponds by convention to the direction of the steady magnetic field in the 
magnet system 2. The measuring co-ordinate system x, y, z to be used may be chosen 
independently of the X, Y, Z co-ordinate system shown in Fig. 1. The gradient coils 3 are 
fed by the power supply unit 4. The MR device also includes an RF transmitter coil 5. The 
25 RF transmitter coil 5 serves to generate RF magnetic fields and is connected to an RF 

transmitter and modulator 6. A receiver coil is used to receive the magnetic resonance signal 
generated by the RF field in the object 7, or in a part of the object 7, to be examined in 
vivo, for example a human or animal body. This coil may be the same coil as the RF 
transmitter coil 5. The magnet system 2 encloses an examination space which is large enough 
30 to accommodate a part of the body 7 to be examined. The RF transmitter coil 5 is arranged 
around or on a part of the body 7 within the examination space. The RF transmitter coil 5 is 
connected, via a transmitter-receiver circuit 9, to a signal amplifier and demodulation unit 
10. The control unit 11 controls the RF transmitter and modulator 6 and the power supply 
unit 4 so as to generate special MR imaging pulse sequences which contain RF pulses and 
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gradients. The phase and amplitude provided by the demodulation unit 10 are applied to a 
processing unit 12. The processing unit 12 processes the applied signal values by way of, for 
example a two-dimensional Fourier transformation, so as to form an MR image. Via an 
image processing unit 13, the MR image is visualized on a monitor 14. 
5 The invention will be described in detail hereinafter, by way of example, 

on the basis of a method which utilizes an MR imaging pulse sequence in the form of a 
known echo planar imaging (EPI) pulse sequence for forming MR images. An EPI pulse 
sequence is known from the cited EP-A-604441. Using these EPI pulse sequences, MR 
signals are measured and an MR image of a part of the body, for example an imaging plane 
10 containing a part of a left ventricle of a heart, is reconstructed from the measured MR 

signals, for example by means of two-dimensional Fourier transformation. However, besides 
said EPI imaging pulse sequences, other known MR imaging pulse sequences can also be 
used, for example a spin echo (SE) imaging pulse sequence, a rapid acquisition and 
relaxation enhancement (RARE) imaging pulse sequence, and a gradient and spin echo 
15 (GRASE) imaging pulse sequence. These imaging pulse sequences are also known from the 
cited EP-A-604441. 

Fig. 2 shows an EPI pulse sequence 20 which includes an excitation RF 
pulse and gradients. The EPI pulse sequence conmiences with the application of an excitation 
RF pulse 100, having a flip angle, and a selection gradient 110 for exciting the spins within a 
20 measuring zone, for example an imaging plane in the x, y plane through the body 7. The flip 
angle amounts to. for example 90°. The selection gradient is oriented in the z direction. 
After the slice selection, an initial phase encoding gradient 120 and a read-out gradient 130 
are applied. Further phase encoding gradients 121, 122, referred to as "blips", are applied 
after the second and further zero crossings of the read-out gradient, so that MR signals 140, 
25 141, 142 can be measured whose sampling instants are simated on lines which are uniformly 
distributed in the k-space. The EPI pulse sequence 30 is repeated for different values of the 
initial phase encoding gradients Gy in order to measure a complete set of MR signals which 
relate, for example 128 or 256 lines in the k-space. In the context of the present patent 
application, a k-space is to be understood to mean a spatial frequency domain in which a 
30 path is traveled during the measurement of the MR signals by application of gradients to the 
steady magnetic field. The path in the k-space is determined by the time integral of the 
gradients applied during the interval between the excitation of the spins until the instant in 
time at which the MR signal is measured. The measured values of the MR signals 
corresponding to the most important pan of the path or paths yield the inverse Fourier 
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transformed values of an image of the imaging plane. 

In order to determine an MR image of an imaging plane which extends 
transversely of a major axis of an ellipsoid approximating the left ventricle of the heart and 
which contains a minor axis of said ellipsoid, the major axis and the minor axis of the 
5 ellipsoid must be estimated. To this end, according to a first version of the method of the 
invention a first MR image of a first imaging plane of the heart and a second MR image of a 
second imaging plane of the heart are formed during a first step, the first and the second 
imaging plane being chosen in such a manner that they extend in parallel and substantially 
perpendicularly to a common axis, for example the head to feet axis, of the body in the 
10 heart, the position of the heart relative to the head to feet axis being estimated in a manner 
customarily used by those skilled in the art. In order to form a first estimate of the major 
axis of the ellipsoid approximating the left ventricle, subsequently a first reference is 
determined from the first MR image and a second reference is determined from the second 
MR image, the first and the second reference corresponding to the same property in the 
15 respective imaging planes of the heart. An example of such a property is a center of a 

circular region in the MR image which represents a transition between the myocardium and 
the blood present in the left ventricle. The major axis of the ellipsoid is subsequently 
estimated on the basis of a connecting line through two points in the steady field which 
correspond to the first reference and the second reference, respectively. The use of said 
20 properly for determining the first estimate of the major axis will be described in detail with 
reference to Fig. 3, 

Fig. 3 shows a first cross-section of the imaging planes and the ellipsoid 
which contains an ellipse 30 and two parallel lines 31 and 32. The ellipse 30 represents a 
cross-section of the left ventricle of the heart in the direction of the major axis. The z axis 
25 represents the major axis of the left ventricle and the t axis represents a minor or 
transverse axis of the left venn-icle. The ellipse 30 is described by the equation 

2 2 

1_ + = 1 where t is a parameter along an axis in the direction of the minor axis, z 
represents a parameter along an axis in the direction of the major axis of the ellipse, and a 
30 and b are two constants which represent the length of the major axis and the minor axis, 

respectively. The parallel lines 31, 32 represent the imaging planes of the left ventricle which 
contain the minor axis of the first and the second MR image, respectively. The parallel lines 
31, 32 are described by the equation z = a^t+^p where is the direction coefficient and p-^ 
is a constant. The points of intersection of a first parallel line 31 and the ellipse 30 are 
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subsequently determined as 



(2) 



A first center 35. having the co-ordinates t^ ^, j and situated between a first and a second 
point of intersection 33, 34 of the first parallel line 31 and the ellipse 30 is determined as 

,0 .„,._f^5^:L m 



b^ . ay 



/?ib 



2 



(4) 



Analogously, a second center 38 is determined between a third and a fourth point of 
15 intersection 36, 37 of the second parallel line 32 and the ellipse 30. A first connecting line 
Ll» constituting the first estimate of the major axis, is given by the equation: 



20 Using the estimate made for the major axis, a third MR image can be formed which 

reproduces a third imaging plane of the left ventricle and is oriented transversely of the first 
estimate LI of the major axis. 

The inaccuracy which occurs in the estimate of the major axis while using 
the first version of the method according to the invention can be classified in two classes. A 

25 first class is governed by the determination of the first connecting line LI and a second class 
is governed by the approximation of the left ventricle by means of the ellipsoid. The 
inaccuracy of the determination of the first connecting line LI follows from the formule (5) 
used and will not be elaborated upon. The inaccuracy of the second class, occurring during 
the estimation of the major axis by the connecting line, can be determined by performing a 
30 first numerical error analysis. According to the first numerical error analysis, the difference 

is determined between the formule (5) for an assumed fixed value of the ratio n = — of 

b ^ 

0.25 and a current value of the formule (5) for a practical value of the ratio n = — which 

a 

may vary from 0.1 to 0.4 and an angulation which may vary from 0"* to 90"*. The result 
of the first numerical error analysis will be described with reference to Fig. 4. 
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Fig. 4 shows a first contour diagram of the inaccuracy as a function of the 
ellipsoid ratio and the angulation. In the present application the ellipsoid ratio is to be under- 
stood to mean the ratio n of the minor axis to the major axis of the ellipsoid and angulation 
is to be understood to mean the angle of the imaging plane of an MR image relative to the 
5 direction of the magnetic field. The surface area below a first line I in Fig. 4 represents a 
first region in which the inaccuracy of the second class remains less than 1 degree in the 
determination of the major axis; the surface area between the first line I and a second line II 
represents a second region in which said inaccuracy remains less than 2 degrees; the surface 
area between the second line 11 and a third line III represents a third region in which said 
10 inaccuracy remains less than 3 degrees; the surface area between the second line III and a 
fourth line IV represents a fourth region in which said inaccuracy remains less than 5 degrees 
and, fmally, the surface area above the line IV represents a fifth region in which said 
inaccuracy amounts to more than 5 degrees. It has been found in practice that said second- 
class inaccuracy is limited to less than I degree for an ellipsoid representing a ventricle of a 
15 hean whose ellipsoid ratio is less than 0.1 with an angulation of the first imaging plane 
amounting to approximately 30 degrees. Said second-class inaccuracy is less than 5 degrees 
for an ellipsoid ratio of 0.25 and an angulation of the first imaging plane amounting to 
approximately 50 degrees. 

In order to make a more accurate second estimate of the major axis, 
20 subsequently a third and a founh MR image are made of a third and a founh imaging plane 
of the heart, respectively, the third and the fourth imaging plane being selected in such a 
manner that they extend parallel to one another and transversely of the first estimate of the 
major axis. Subsequently, a third and a fourth reference are determined from the third and 
the fourth MR image, respectively, said references corresponding to a third and a fourth 
25 center, respectively, of circular regions in the third and the fourth MR image, respeciively, 
which represent the transitions between the myocardium and the blood. A second estimate of 
the major axis is then determined by a second connecting line L2 through two points in the 
steady magnetic field which correspond to the third and the fourth reference, respectively, in 
the same way as described with reference to Fig. 3. 
30 Using, for example the first version of the method according to the 

invention, in practice a first estimate of the major axis can be made on the basis of an 
angulation of the fu-st and the second MR imaging plane which amounts to 50 degrees and an 
ellipsoid ratio of 0.1, an inaccuracy of the second class amounting to approximately 5 
degrees for the fu-st estimate of the major axis. Subsequently, an inaccuracy of the second 
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estimate of the major axis of approximately 1 degree can be achieved by way of a second 
estimate of the major axis which is determined by repeating the steps of the first version for 
the third and the fourth MR image. This method can be used in practice for the estimation of 
the major axis of an ellipsoid which represents the left ventricle of the hean with an ellipsoid 
5 ratio which may vary in a range from 0.1 to 0.4 and an angulation of the imaging planes 
which may vary in a range from 0 to 70 degrees. 



of the ellipsoid, use is made of the three-dimensional geometrical properties of an ellipsoid of 
revolution so as to determine also the orientation of the minor axis relative to the steady 
10 magnetic field. This version of the method will be described in detail with reference to Fig. 
5. 



ellipsoid, which second cross-section contains the ellipse 30 and the parallel lines 31 and 32. 
Fig. 5 also shows the coimecting line LI and an orientation line M. The first connecting line 

15 LI between the first center 35 and a second center 38 of the first and the second point of 
intersection 33, 34, respectively, of the first parallel line 31 with the ellipse 30 and the third 
and the fourth point of intersection 36, 37, respectively, of the second parallel 32 is 
determined by means of the formules (1), (2), (3), (4) and (5). The orientation line M is 
directed substantially perpendicularly to the parallel lines 31, 32. The angle between the 

20 orientation line M and the minor axis of the ellipsoid t is then given by the equation: 



in which ^ represents the angle between the direction of the orientation line M and the minor 
25 axis of the ellipsoid t, and represents the angle between the first connecting line LI and the 
orientation line M. The direction of the orientation line M is known because it has been 
chosen to extend perpendicularly to the first and the second imaging plane of the left 
ventricle and the directions of the first and the second imaging plane relative to the direction 
of the steady magnetic field are known. The direction of the first connecting line LI is 
30 known from the first estimate of the major axis of the ellipsoid. An improved estimate of the 
major axis is offered by rotation of the orientation line M in the direction of the first 
coimecting line LI through an angle amounting to 7r/2-f. The validity of the formule (6) is 
limited to a positive value of a discriminant wherefrom a square root is determined. In the 
present example this means that the maximum angle c^^, between the first connecting line LI 



For further enhancement of the accuracy of an estimate of the major axis 



Fig. 5 shows a second cross-section of the imaging planes and the 



f = tan 



-1 



2 



(6) 
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and the orientation line M, corresponding to the maximum angle between planes of the 
imaging planes of the heart to be imaged and a major axis of the ellipsoid which represents 
the left ventricle of the heart of a body in the examination space, is determined by the 
formule 




(7) 



For an ellipsoid ratio n for a human body, customarily having a value of, 
for example 0.25, formule (7) then yields a value of 60 degrees for the iimu angle for the 
10 angulation <^c* 

An inaccuracy of the first class in the second estimate of the major axis 
while using the second version of the method of the invention is then determined by the 
formule (6) and the absolute inaccuracy in S and will be explained with reference to Fig. 6. 
Fig. 6 shows the effect of an inaccuracy of the first class in a line 60 which represents the 

15 function for the range 0 ^ 4> ^ where 0, = 1.08L The angle c^> is plotted in 

radians along the horizontal axis and the angle f is plotted in radians along the venical axis. 
Fig. 6 shows that an inaccuracy in the determination of the angle 0 yields a substantially 
equal inaccuracy in the determination of the angle The inaccuracy of the second class, 
occurring during the estimation of the major axis by the connecting line while using the 

20 second version, can be determined by performing a second numerical error analysis. 

According to the second numerical error analysis the difference is determined between the 

formule (6) for an assumed fixed value of the ellipsoid ratio n = _ amounting to 0.25 

^ b 

and a current value of the formule (6) for a practical value of the ellipsoid ratio n = _ 
which may vary from 0.1 to 0.4 and an angulation a, which may vary from 0 to 90 degrees. 
25 The result of the second nimierical error analysis will be described with reference to Fig. 7. 

Fig. 7 shows a second contour diagram of the inaccuracy as a function of 
the ellipsoid ratio n and the angulation. The second contour diagram shows the lines 1, II, II, 
IV which define the first, the second, the third, the fourth and the fifth region. The value of 
the inaccuracy in said regions is the same as described wi± reference to Fig. 5. 
30 It is also possible to determine the first estimate of the major axis by 

means of the first version of the method and to determine subsequently the third and the 
fourth MR image, followed by the determination of a more accurate second estimate of the 
major axis by means of the second version of the method. 

A known edge detection algorithm can be used to determine a center of a 
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circular region in an MR image of an imaging plane of the left ventricle. This algorithm is 
known from the anicle "Discrete Dynamic Contour Model", published by S. Lobregt and 
M.A. Viergever in IEEE Transactions on Medical Imaging, Vol. 14, No. 1 (1995), pp. 12- 
24. Subsequently, a center of the circular region is determined by way of a center detection 
5 algorithm, said center representing the transition between the myocardium and the blood in 
the left ventricle of the heart. 

In order to enhance the accuracy of the edge detection algorithm, the 
contrast between the part of the MR image representing the myocardium and the pan 
representing the blood in the left ventricle can be increased. To this end, for example use can 
10 be made of a TjY pre-pulse sequence. The T^Y pre-pulse sequence is known from the cited 
anicle "Myocardial Suppression In Vivo by Spin Locking with Composite pulses" , published 
by W.T. Dixon et al. The T^Y pre-pulse sequence is generated prior to the EPI pulse 20. 
Said contrast amounts to approximately 2 in an MR image formed in a customary manner by 
means of, for example EPI pulse sequences. When a TjY pre-pulse sequence is applied prior 
15 to the EPI pulse sequences, said contrast in the MR image is increased to 4. 

Another possibility for a pre-pulse sequence preceding the EPI imaging 
pulse sequence is the generating of a known To pre-pulse sequence prior lo the EPI pulse 
sequence 20. The Tj pre-pulse sequence is known from the article "Coronary Angiography 
with Magnetization- Prepared T2 Contrast", published by J.H. Brittain et al. in Magnetic 
20 Resonance in Medicine, No. 33, pp. 689-696, 1995. The T2 pre-pulse sequence is generated 
prior to the generating of the excitation RF pulse 100 of the EPI pulse sequence 20 and 
suppresses muscular tissue and vascular tissue in an MR image of, for example the left 
ventricle of the heart. 

Another possibility for enhancing the contrast is the introduction of a 
25 contrast medium into the blood stream, so that in an MR image of an imaging plane of the 
heart the blood contained in the myocardium is reproduced with a higher contrast than the 
myocardium. Gd-DPTA is an example of a known contrast medium. For such a contrast 
medium the measured T^ and T2 are influenced in conformity with 



30 



t;^ = ^2' - ^2'^ 

where Rj, Ro are respective constants and C is the concentration. For a magnetic field 
strength of, for example 1.5 T, the values of die two constants R^, Rj are approximately 4.5 
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ms"^ (minol/kg)'^ 
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CLAIMS: 



1. 



A magnetic resonance (MR) imaging for imaging a pan of a human or 



animal body arranged in a steady magnetic field, which method involves the use of imaging 
pulse sequences for the measurement of the MR signals for reconstructing the MR images, 
approximation of the part of the body by an ellipsoid, imaging of an imaging plane of the 
5 part which contains a minor axis of the ellipsoid, which minor axis is derived from a major 
axis of the ellipsoid, and 

which method includes the following steps for estimating the major axis: 

a) formation of a first MR image of a first imaging plane and a second MR image 
of a second imaging plane of the part of the body, the imaging planes extending transversely 

10 of a coTTunon axis of the body. 

b) determination of a first reference from the first MR image and a second 
reference from the second MR image, and 

c) determination of a first estimate of the major axis of the ellipsoid from a 
connecting line through two points which correspond to the first reference and the second 

15 reference, respectively, 

characterized in that the first and the second reference are determined by one and the same 
property. 

2. A method as claimed in Claim 1, characterized in that the method also 

20 includes steps for 



of a fourth imaging plane of the pan of the body, the third and the founh imaging plane 
being chosen so as to extend transversely of the first estimate made for the major axis, and a 
step in which a second estimate is made for the major axis by repeating the steps b) and c) 
25 while replacing the first MR image by the third MR image and the second MR image by the 
fourth MR image. 

3. A method as claimed in Claim 1 or 2, characterized in that said property 

is determined by a center in a circular region of said MR images of the imaging planes of the 



formation of a third MR image of a third imaging plane and a fourth MR image 
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part of the body. 

4 A method as claimed in Claim 1, 2 or 3, characterized in that in order to 

make a next estimate of the major axis which replaces a previous estimate of the major axis 
5 the method also includes the following steps 

determination of an orientation line through a first point which corresponds to 
the first reference, substantially perpendicularly to the imaging planes, and 

determination of the next estimate of the major axis from a combination of the 
predetermined ratio of the minor axis to the major axis of the ellipsoid and an angle between 
10 the previous estimate of the major axis and the orientation line. 

5^ A method as claimed in Claim 4, characterized in that the predetermined 

ratio of the minor axis to the major axis of the ellipsoid is chosen to be approximately 1/4. 

15 6. A method as claimed in Claim 1, 2, 3, 4 or 5, characterized in that the 

imaging pulse sequences include an Echo Planar Imaging pulse sequence. 

7. A method as claimed in Claim 1, 2, 3. 4, 5 or 6, characterized in that in 
order to enhance the contrast between the parts representing blood and the parts representing 

20 tissue in said MR images, the method includes a step in which a pre-pulse sequence is 
generated prior to the generating of the MR imaging pulse sequence. 

8. A method as claimed in Claim 1.2. 3, 4, 5, 6 or 7, characterized in that 
in order to enhance the contrast in said MR images the method includes a step in which a 

25 contrast medium is introduced into a blood stream through the pan to be imaged. 

9. A magnetic resonance (MR) device for imaging a part of a body of a 
human or animal which is to be arranged in the device, which MR device also includes: 
means for sustaining the steady magnetic field. 

30 means for generating RF pulses, 

means for generating gradients in the steady magnetic field, 

a control unit for generating control signals for the means for generating RF pulses and the 
means for generating the gradients, 
means for measuring MR signals. 
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means for reconstructing a first MR image and a second MR image from the measured MR 
signals from a first imaging plane and a second imaging plane, respectively, of a part of the 
body to be arranged in the device, which imaging planes are directed transversely of a 
conmion axis of the body to be arranged in the device, 
5 means for determining a first reference from the first MR image of the first imaging plane, 
means for determining a second reference from the second MR image of the second imaging 
plane, 

me£ins for determining a first estimate of a major axis of an ellipsoid which approximates the 
part of the body to be arranged in the device, said first estimate being determined by a 
10 connecting line between two points, said two points corresponding to the first reference and 
the second reference, respectively, and 

means for forming an MR image of an imaging plane of the pan which contains a minor axis 
of the ellipsoid, 

characterized in that the means for determining the first reference and the means for 
15 determining the second reference are arranged to determine the same property. 
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FIG.1 



BNSOOCfO: <WO 991S450A1_I_> 



wo 99/18450 



PCT/IB98/01496 




BNSOOaO; <V/0 99ie450Al_1_> 




BNSCOCID: <WO_9918<50A1J_> 



INTUUNA riONAK SKARCIl RKPOR I 



hUcnintiuiial application No. 

PCT/IB 98/01496 



A. CLASSU'ICA TION OK SUIUHCI MA TI I:K 



IPC6: GOIR 33/56, A61B 5/055 ^ ^ . 

According lo International Patenl Ctoifiealion (IPC) «r Ui hoih naUona? cUs.^irn:al»in and II C 



Minimum documentation searched (clas.tiricaiiiin system followed by classincalion symbols) 

1PC6: GOIR, A618 



Documentalion searched other than minimum dncumcnUOnn to Oic extent that .such documcnus are included in the neld.s searched 

SE,DK,FI,NO classes as above 



Electronic data base consulted during the international search (name uf data bhse and, where pracltcahk, .<tcarch terms used) 



WPI, PA J, INSPEC 



C DOCUMUKIS CONSIDIiRI'D'I O BH RI-IJ-VAN J 



Category* 



Citation of di>cunient, with indicntton, where approprinlc, of the i clcvaut pas.sngcs 



US 5680471 A (TOYGMITSU KANEBAKO ET AL) , 
21 October 1997 (21.10.97), column 1, 
line 16 - column 2, line 59, figures 4,5, 
abstract 



US 5107838 A (KOJIRO YAMAGUCHI) , 28 April 1992 
(28.04,92), column 2, line 5 - line 29 



US 5427100 A (NOBORU HIGASHI ET AL) , 27 June 1995 
(27.06.95), figures lA-lC, abstract 



Relevant to cinini No. 



1-9 



1-9 



1-9 



[ Turlher documents arc Hslcd in lite coiitinjintioii ol" Hox C. [ xl p.it*-Mii fainily annex. 



♦ Special cateijorics of cited documcnl-T 

"A' document defining the ^cral stale of the art \vhich i% noi cfinsidered 

lo he of particular relevance 
"fi" criicr document but published on or after the tnlcmalinnal filing dale 

'L* document which may throw duuht^ on pritirity claim(s) or ^vhich is 
cited Ui establish the publication date of antjlhcr cilatinn or other 
special reason (as jpcdfied) 

"O' document referring to an nral disclosure, tt?e, cxhihiiion or olhcr 
means 

*P' document published prior to ihc intemaliona] filing dale but later than 
the priority dale d aimed 



""I"* IMcr dncumcnl puhlijrticd atlcr Vhc micmaiional filing dale or pricinly 
d;»lc and not in crmJIicl wlh the ftppiic*tion bul alcd \o undcrsUind 
the principle nr thctjry undcrl>ing the invention 

*X" d<K*umenl of particular relevance: the claimed mvcniion cannoi he 
ctmsidcrcd novd or canntrt b« considered l»» involve an inventive 
sicp when the documem is iRltcn altmc 

' Y" dticument of particular rclc\-ancc: the claimed invcnlion cannot be 
con'adcrcd to involve an inventive step when the document is 
combined with one or more other such documents, such wmhinaiion 
being obvious lo a person skilled in lb e art 
'&" diicumtnt member of the same patent family 



Dale of mailing of iht: iiiicriiaiional search report 



Dale ol thc actual coinplcUon of Uic international search 



IS March 1999 



Name and mailing address of Ihc ISA/ 

Swedish Patent Office 

Box 5055, S-102 42 STOCKHOLM 

Facsimile No. + 46 8 66f» 02 RC 



1 6 -03- 1999 



Anthorir.cd oHiccr 

Joni Sayeler 

•relep1t"i*c No. + 46 8 7K2 25 00 



Form PCl7ISA/2tO (second sheet) (July 1992) 

BNSCXXID: <V/0 991B450A1J_> 



2 

IM KRNA nONAK SI- ARCH HKPOR I 



Inlcnintional application No, 

PCT/IB 98/01496 



C (Conlinuntimi). IJOCUM liN I S C;ONSIl)i:Hi:n I () HI- Rin.IA'AN I 



Category^ 


Citalion of di)CUfuciiV. willi iiiUicnliim, w1>l'i c npprnpriale, of llic rdcvnnl pnssngcs 


RclevaiU to clnini No. 


A 


US 5360006 A (EDWARD A. GEISER ET AL) . 

1 November 1994 (01.11.94). column 3, 
line 32 - column 4, line 32 


1-9 



Form PCr/ISA/210 (conlinualion of second sheet) (July I9»J2) 

DNSEXX;it>. <WO 9918+50A1 J_> 



t 



IN TKRNA I IONAL SKARCII RICPOU 1 

InforiTialioii on patent fniiiily members 02/03/99 


[iilei iiaiitMial applic*Tiio!i Ko. 

PCT/IB 98/01496 


I'aieni dticumcm 
ciled in search rep on 


I'ubiicaUon 
dale 


Palcm familv 1 Publivalion 
membcr(sj | ^a^^* 



US 


5107838 


A 


28/04/92 


JP 


4218139 


A 


07/08/92 


US 


5427100 


A 


27/06/95 


JP 


5212013 


A 


24/08/93 


US 


5360006 


A 


01/11/94 


EP 
JP 
WO 


0533782 
6503483 
9119457 


A 
T 
A 


31/03/93 
21/04/94 
26/12/91 



Form PCt71SA/2lO (palcnl family annex) (July IVV2) 

BNSDOCID: <WO 991B45QAl_l_> 



